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FINITE-ELEMENT  ANALYSIS  OF  LAMINATED 
BIMODULUS  COMPOSITE-MATERIAL  PLATES 

J.N.  REDDY  and  W.C.  CHAO 

School  of  Aerospace,  Mechanical  and  Nuclear  Engineering 
The  University  of  Oklahoma,  Norman,  OK  73019 

Abstract  - Fi ni te  element  analysis  of  the  equations  governing  the  small- 
deflection  elastic  behavior  of  thin  plates  laminated  of  anisotropic  bi- 
modulus materials  (which  have  different  elastic  stiffnesses  depending 
upon  the  sign  of  the  fiber-direction  strain)  is  presented.  Single-layer 
and  two-layer  cross-ply,  simply-supported  rectangular  plates  subjected  to 
sinusoidally  distributed  normal  pressure  and  uniformly  distributed  normal 
pressure  are  analyzed.  For  the  sinusoidally  loaded  case,  the  finite  ele- 
ment solutions  are  gratifyingly  close  to  the  exact  closed-form  solutions. 

NOMENCLATURE 

A,-  = stretching  stiffness 

1 J 

a,b  = dimensions  of  plate  in  x and  y directions 
B- . = bending-stretching  stiffness 

c = a/b 

D^j  = bending  stiffness 

c t 

En»En>etc.  = respective  compressive  and  tensile  Young's  moduli  for 
orthotropic  bimodulus  material  (without  subscripts  for 
isotropic  bimodulus  material) 

F.j  = finite-element  force  components 
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shear  moduli  for  orthotropic  bimodulus  material 

total  thickness  of  plate 

E22Wh1 * 3/qa4 

stiffness  coefficients  in  the  finite-element  formulation 
stress  couples  and  stress  resultants  (i =1,2,6) 
plane-stress-reduced  stiffnesses  for  orthotropic  material 
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= -normal  pressure 
= maximum  normal  pressure 
= midsurface  of  the  plate 
= finite-element  matrix  coefficients 
= displacements  in  x,y,z  directions 

= plate  coordinates  in  longitudinal,  transverse,  and  downward 
thickness  directions 
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‘ !nx/h- V 

= neutral -surface  positions  associated  with  =0  and  e = 0 

x y 

= penalty  parameters 
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= fiber- direction  strain  in  z-th  layer 

= strain  component  at  arbitrary  location  (x,y ,z) ,( j = l ,2,6) 

= finite-element  interpolation  functions 
« stress  component  ( j -1  ,2,6) 

= principal  Poisson's  ratios  for  orthotropic  bimodulus  material 

* angle  between  fiber  direction  and  plate  reference  direction 
(x  axis) 

= slopes,  3w/3x , 3w/3y  (and  slope  functions  in  the  YNS  theory) 


1 . INTRODUCTION 

Fiber-reinforced  composite  materials  are  finding  widespread  use  in 

many  engineering  structures.  This  is  primarily  due  to  their  high  stiffness- 
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to-weight  ratio,  and  anisotropic  material  properties  that  can  be  tailored 
through  variation  of  the  fiber  orientation  and  stacking  sequence.  Material 
characterization  of  composite  materials  has  shown  that  certain  composite 
materials,  known  as  bimodulus  materials,  behave  differently  in  simple  tension 
and  compression  under  static  loading  [1].  The  analysis  of  bimodulus  mater- 
ials is  more  difficult  than  that  of  ordinary  materials  (i.e.,  materials  whose 
elastic  moduli  are  the  same  in  tension  and  compression)  since  the  elastic 
moduli  depend  on  the  sign  of  the  stress,  which  is  unknown  a-priori.  A plate 
bent  by  transverse  load  experiences  both  tension  and  compression,  and  there- 
fore the  analysis  of  plates  made  of  bimodulus  material  becomes  more  difficult. 

The  first  material  model  appropriate  for  bimodulus  materials  is  appar- 
ently due  to  Ambartsumyan  and  Khachatryan  [2,3],  There  have  been  only  a small 
number  of  bimodulus  material  applications  to  structural  plate  problems  (see 
[4-11]).  All  of  these  analyses  were  limited  to  isotropic  materials  with 
different  properties  in  tension  and  compression.  Ambartsumyan  [3]  and  Tabaddor 
[9]  proposed  linearized  material  models  for  bimodulus  anisotropic  materials, 
whili  Jones  and  Morgan  [10]  presented  a closed-form  solution  for  cylindrical 
bending  of  a thin,  cross-ply  laminate.  More  recently,  Bert  [11]  presented 
a macroscopic  material  model  for  bimodulus  fiber-reinforced  composites.  This 
model  was  found  to  agree  well  with  experimental  results.  A summary  of  various 
bimodulus  material  models  can  be  found  in  the  survey  paper  of  Bert  [12]. 

Using  the  material  model  of  Bert  [11],  closed-form  solutions  have  been  devel- 
oped for  clamped  elliptic  plates  under  uniform  loading  and  simply-supported 
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The  present  paper  is  concerned  with  the  finite-element  analysis  of 
thin  elastic  plates  laminated  of  anisotropic  bimodulus  materials  described 
by  the  model  of  [11].  Numerical  results  are  presented  for  single-layer 
and  two-layer  cross-ply,  simply-supported  rectangular  plates  under  both 
sinusoidal  loading  and  uniform  loading.  Only  in  the  case  of  sinusoidal 
loading  does  a closed-form  solution  exist  [13]  for  rectangular  plates.  In 
this  case, the  present  finite-element  solution  closely  agrees  with  the  exact 
closed-form  solution.  Apparently,  the  present  investigation  is  the  first 
to  consider  the  finite-element  analysis  of  anisotropic  bimodulus  plates 
of  finite  dimensions. 


2.  THEORY  AND  FORMULATION 

Consider  a plate  of  constant  thickness  h composed  of  thin  anisotropic 

layers  oriented,  generally  speaking,  at  angles  0 1 , 0 2 The  origin  of  the 

coordinate  system  is  located  within  the  middle  plane  (x,y)  with  the  z-axis 
being  normal  to  the  midplane.  The  material  of  each  layer  is  assumed  to 
possess  a plane  of  elastic  symmetry  parallel  to  the  xy-plane.  The  thin- 
plate  equilibrium  equations,  in  the  absence  of  body  forces  and  body  moments, 
are 


3Ni  3N6  3N6  3N2 


All  of  the  symbols  have  the  usual  meaning  and  are  defined  in  the  Nomenclature. 

In  the  thin-plate  theory  (Kirchhoff  hypothesis),  the  strain-displace- 
ment relations  are  given  by 
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where  0 and  0 are  the  slopes, 
x y 


9 - |W  = o f 9 - f*=  0 
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(2.4) 


Equations  (2.1 )-(2.4)  are  valid  for  both  ordinary  and  bimodulus  mat- 
erials. What  distinguishes  a bimodulus  material  from  an  ordinary  material 
are  the  constitutive  relations.  Following  Bert's  [11]  fiber-governed  macro- 
scopic material  model,  we  assume  that  there  are  two  symmetric  plane-stress 
reduced  stiffness  matrices:  one  when  the  fibers  are  in  tension  along  their 
length  and  another  when  they  are  in  compression  in  the  same  direction.  Then 
the  stress-strain  relation  for  a thin  orthotropic  bimodulus  material  may  be 
written  as 
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where  Q..^  denote  the  plane-stress-reduced  stiffnesses  for  2.-th  aniso- 
tropic layer  in  tension  (k=l)  or  compression  (k=2): 
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where  denotes  the  fiber-direction  strain  at  any  arbitrary  point  in  layer, 


The  stretching,  bending-stretching  coupling,  and  bending  stiffnesses  of 
the  bimodulus  laminate  are  defined  exactly  in  the  same  way  as  in  the  case  of  an 
ordinary  laminate. 
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<WV  ' <]"  Q1  jti0 -z-z2  )dz  + I W C.z.z^dz)  (2.7) 


where  z is  the  distance  from  the  midsurface  to  the  neutral  plane  (which  is 
n 

unknown  a-priori ) . Note  that,  in  addition  to  performing  the  integrations 

in  a piecewise  manner  from  layer  to  layer,  one  must  take  into  account  the 

possibility  of  different  properties  (tension  or  compression)  within  a layer. 

For  example,  in  the  two-layer,  cross-ply  case,  the  coefficients  A.,  are 

■ J 

given  by 
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Q(Jn  dz  (2.8) 


wherein  it  is  assumed  that  layer,  2 = 1 (bottom)  occupies  the  thickness 
space  from  z = 0 to  z * h/2  and  layer,  2=2  (top)  occupies  the  thick- 
ness space  from  z = -h/2  to  z = 0 (z  is  measured  positive  downward  from  the 

midplane).  In  writing  eqns.  (2.8)  it  is  assumed  that  the  upper  portion  of 
the  top  layer  (2=2)  is  in  compression  (k=2)  in  the  fiber  direction  and  that 
the  lower  portion  of  the  top  layer  is  in  tension  (k=l),  while  the  portion 
z = 0 to  z = znx  of  the  bottom  layer  (2=1)  is  in  compression  and  the  por- 
tion z = z to  z = h/2  is  in  tension.  Here  z and  z denote  neutral- 

nx  nx  ny 

surface  positions  associated  with  e = 0 and  e = 0,  respecti vely. 

x y 
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Now  the  laminate  constitutive  relations  for  an  arbitrary  laminate  of 
anisotropic  material  are  (using  eqns.  (2.2),  and  (2.5)), 
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Substituting  eqns.  (2.9)  into  eqns.  (2.1),  one  obtains  the  governing' 
equations  in  terms  of  the  displacements.  For  the  finite-element  formulation, 
we  construct  the  associated  functional.  The  problem  of  seeking  solutions 
to  eqns.  (2.7)  and  (2.9)  is  equivalent  to  minimizing  the  total  potential 
energy, 

rh/2 
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subject  to  the  contraints  in  eqn.  (2.4).  Note  that  if  9 and  9 are 

x y 


eliminated  from  eqn.  (2.10)  using  relations  (2.4),  the  resulting  total 
potential  energy  functional  involves  the  second-order  derivatives  of  the 
transverse  deflection,  and  consequently  imposes  severe  continuity  require- 
ments on  the  finite-element  interpolation  functions. 

To  include  the  constraint  conditions  into  the  variational  formulation, 
we  use  the  penal  ty -function  method  (see  Reddy  [14]).  In  the  penalty  method, 
a constrained  minimization  problem  is  replaced  by  an  unconstrained  minimi- 
zation problem  whose  solution  converges  to  the  true  solution  in  the  limit 
certain  parameter,  called  the  penalty  parameter,  approaches  the  value  of 
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infinity.  For  details,  the  reader  is  referred  to  the  senior  author's  recent  paper 
[14]. 

The  modified  functional  to  be  minimized  is  given  by 


ifp(  U , V ,W,9^  ,9^ ) = rr(u,V,W,9x,9y)  + 


[ai  (ll-  + a2  9J^  dxdy  (2.11) 
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where  and  a2  are  the  penalty  parameters  whose  values  are  to  be  preassigned 
(of  the  order  1012  - 1015).  Generally  speaking,  the  penalty  functional  tt 


approaches  tt  as  aj  and  a2  approach  «,  and  for  finite  values  of  and  ct2, -r 


has  no  physical  significance.  However,  in  the  case  of  plates,  tt  has  a 


physical  significance  (only  incidental):  for  the  following  values  of  and 
a2,  it  corresponds  to  the  functional  associated  with  the  Yang,  Norris,  and 
Stavsky  theory  (a  generalization  of  Mindlin's  theory  for  isotropic  plates  to 
laminated  anisotropic  plates)  [14]: 
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where  ki  and  k2  are  the  shear  correction  factors,  and  9 and  9 are  the 

* y 


slope  functions  of  the  thick-plate  theory  (not  equal  to  the  slopes  of  w). 

Thus  (in  retrospect  only)  the  Mind! in  thick-plate  theory  can  be  interpreted 
as  a theory  obtained  from  the  thin-plate  theory  by  perturbation;  for  suffi- 
ciently large  a's  one  recovers  the  thin-plate  theory  from  the  thick-plate 
theory. 


3.  FINITE-ELEMENT  MODELS 


We  now  construct  a finite-element  model  based  on  t (u,v,w,9  ,9  ). 

P x y 


We  assume,  over  each  element  R , the  same  kind  of  interpolation  for  all  of 


the  variables, 
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S ^ p p 0 ^00 

U = Z Uj  , v = z vi  , etc.  (n  = nodes  per  element)  (3.1) 

where  are  the  element  interpolation  (or  shape)  functions,  and  u®,  and 
v®  are  the  nodal  values  of  ue  and  ve,  respectively.  Substituting  (3.1) 
into  the  first  variation  of  t (u,v,w,9  ,9  ),  and  collecting  the  coeffi- 

p x y 

cients  of  the  variations,  5u^ , <Sv.,  etc.,  we  obtain 
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where  Ue}  = { { u e } { ve } {we}  C 8 e > {9®}}',  and  the  elements  K^6  (a, 3=1, 2, 
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...,5;  i ,j=l  ,2,. . . ,n)  of  the  stiffness  matrix  are  given  by 
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The  element  stiffness  matrices  are  assembled  in  the  usual  manner,  and 
boundary  conditions  of  the  problem  are  imposed  before  solving  for  {A}. 

In  the  present  study  linear  (n=4)  elements  of  the  serendepity  family 
are  used.  The  element  stiffness  matrices  for  the  linear  element  are  of 
order  20x20. 

4.  NUMERICAL  RESULTS 

In  the  following  we  present  numerical  results  for  (a)  single-layer, 
isotropic  and  orthotropic  bimodulus  rectangular  plates,  and  (b)  two-layer, 
cross-ply  (0°/90° ),  orthotropic  bimodulus  rectangular  plates  both  with  simply- 
supported  edge  conditions.  A mesh  of  3x5  in  the  quarter  plate  is  employed. 

The  side-to- thickness  ratio  is  taken  to  be  100.  For  simply-supported  plates 
under  sinusoidal  loading,  the  neutral  surfaces  are  planes  and  closed-form 
solutions  exist  [13]. 

For  single-layer,  isotropic  bimodulus  rectangular  plate,  approximate 
solution  (Kamiya  [7]),  and  closed-form  solutions  (Bert  [13])  are  available. 
Table  4.1  presents  a comparison  of  dimensionless  maximum  deflections  obtained 
by  various  methods  for  the  case  of  sinusoidal  loading.  Dimensionless  neutral- 
surface  position  for  the  same  case  are  compared  in  Table  4.2.  Numerical 
results  obtained  by  all  of  the  methods  are  very  close.  Note  that  the  loca- 
tion of  the  neutral  surface  is  independent  of  the  plate  aspect  ratio. 
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Table  4.1.  Comparison  of  dimensionless  maximum  deflections  for 
sinusoidally  loaded,  simply  supported  square  plate 
made  of  isotropic  bimodulus  materials  (vc=0.20,  b/h=l 00) 


Finite-element  solution 
(present  work) 


1 0.5 

! i.o 

i 2.0 


— 

0.5 

0.04387 

1.0 

1 .0 

0.02956 

2.0 

0.02040 

0.5 

0.11230 

2.0 

1.0 

0.07569 

2.0 

0.05223 

Approximate 
solution  [7] 


« = wmax  (x=0’^0) 


0.007018 

0.004730 

0.003264 


0.0435 


0.0194 


0.04351 

0.02933 

0.02024 

I 


! 0.11139 
| 0.07508 
! 0.05182 


0.2833 

0.0 

0.3224 


0.1417 

0.0 

0.1928 


0.2833 

0.0 

0.3224 


0.2835 

0.0 

0.3226 


Table  4.2.  Values  of  the  dimensionless  neutral -surface  position  for  the 
isotropic  bimodulus  rectangular  plate  of  Table  4.1  ( b/h=l 00 ) 


Et/Ec 


a/b  = 0 

.5 

a/b  = 

.0 

a/b  = 

2.0 

Closed-form 
sol ution[13] 

Present 

FEM 

sol uti on 

Closed-form 
sol ution[13 ] 

Present 

FEM 

solution 

Closed-form 
sol ution[13 ] 

Present 

FEM 

solution 

- 0.08951 

0.0 

0.10138 

1 

- 0.08951 

0.0 

0.10189 

-0.08951 

0.0 

0.10188 

- 0.08951 

0.0 

0.10189 

-0.08951 

0.0 

0.10188 

- 0.08951 

0.0 

0.10189! 
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For  the  orthotropic  bimodulus  plates,  the  following  elastic  properties 


(see  Bert  [13])  are  used: 


Arami d-Rubber 

Tensile 

Properties 

Compress 

Properti 

Major  Young's 

modulus 

(GPa) 

3.5842 

0.0120 

Minor  Young's 

modulus 

(GPa) 

0.00909 

0.0120 

Major  Poisson' 

's  ratio 

0.416 

0.205 

Shear  modulus 

(GPa) 

0.0037 

0.0037 

Polyester-Rubber 

Major  Young's 

modul us 

(GPa) 

0.617 

0.0369 

Minor  Young's 

modulus 

(GPa) 

0.008 

0.0106 

Major  Poisson 

's  ratio 

0.475 

0.185 

Shear  modulus 

(GPa) 

0.00262 

0.00267 

Table  4.3  shows  the  closed-form  solutions  [13],  and  the  finite-element 
solutions  for  dimensionless  deflections  and  neutral -surface  locations  of 
single-layer,  simply-supported  rectangular  plate  subjected  to  sinusoidal 
loading.  There  is  excellent  agreement  between  the  closed-form  solutions 
and  the  finite-element  solutions. 

Numerical  results  are  presented  in  Table  4.4  for  two-layer,  cross-ply 
laminated  plate  (0°/90°)  subjected  to  sinusoidal  loading.  Again,  the  closed- 
form  results  and  the  finite-element  results  are  gratifyingly  close  for  both 
the  deflections  and  neutral -surface  locations  (which  are  independent  of  x 
and  v ) . 

The  example  problems  considered  thus  far  admit  closed-form  solutions 
due  to  the  fact  that  the  neutral-surface  locations  are  constant  (i.e., 
independent  of  x and  y)  for  the  particular  geometry,  support  conditions,  and 
loading  considered.  For  the  same  geometry  (i.e.,  rectangular)  and  boundary 
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Polyester-Rubber  Aramid-Rubber  Material 
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Table  4.3.  Fiber-direction  neutral-surface  location  and  deflections 

for  simply-supported  rectangular  plate  of  single-layer  (0°) 
aramid-rubber  and  polyester-rubber  (b/h  = 100) 


Closed- form 


0.4457 

0.4457 

0.4444 


0. 

0. 


0.4368 

0.4334 

0.4298 


0.3040 

0.3040 

0.3040 

0.3040 

0.3040 

0.3027 

0.3014 

0.2998 

0.2979 

0.2958 

0.2936 


nx/h 

w = w Ec 

max  22 

hVq0b4 

Finite-element  solution 

Fi  ni te- 

Closed-form 

Fi ni te- 

element 

solution 

solution[13] 

element 

solution 

u=103u  /w 

max  max 

v=103v  /w 

max  max 

0.4454 

0.001881 

0.001875 

0.2808 

0.1306 

0.4451 

0.003661 

0.003640 

0.2342 

0.1093 

0.4447 

0.006253 

0.006211 

0.2007 

0.09171 

0.4440 

0.009679 

0.009605 

0.1754 

0.07738 

0.4431 

0.01387 

0.01376 

0.1556 

0.06583 

0.4420 

0.01870 

0.01854 

0.1397 

0.05647 

0.4393 

0.02956 

0.02928 

0.1158 

0.04270 

0.4363 

0.04089 

0.04049 

0.09864 

0.03337 

0.4328 

0.05170 

0.05120 

0.08562 

0.02686 

0.4292 

0.06143 

0.06085 

0.07548- 

0.02222 

0.4254 

0.06995 

0.06931 

0.06732 

0.0T882 

0.3041 

0.000816 

0.000821 

0.1974 

0.08969 

0.3041 

0.001655 

0.001656 

0.1588 

0.07239 

0.3039 

0.002975 

0.002968 

0.1364 

0.06110 

0.3036 

0.004888 

0.004866 

0.1194 

0.05206 

0.3031 

0.007478 

0.007434 

0.1061 

0.04482 

0.3026 

0.01079 

0.01071 

0.09543 

0.03966 

0.3012 

0.01954 

0.01935 

0.07930 

0.03062 

0.2995 

0.03056 

0.03021 

0.06764 

0.02495 

0.2976 

0.04278 

0.04224 

0.05882 

0.02102 

0.2954 

0.05505 

0.05431 

0.05192 

0.01821 

0.2931 

0.06652 

0.06560 

0.04636 

0.01616 

Polyester-Rubber  Aramid-Rubber  I Material 
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Table  4.4.  Neutral-surface  positions  and  deflections  for  simply- 

supported  rectangular  plates  of  cross-ply  laminated  aramid- 
rubber  and  polyester-rubber  (sinusoidal  loading  case). 


Aspect 
ratio  c 

Z 

X 

Ll 

WE22  h3/qb4 

Closed-form 
sol uti on[l 3] 

Finite- 

element 

solution 

Closed- form 
sol  uti on[l 3] 

Fini te- 
element 
sol uti on 

Closed-form 

solution[13] 

Fini te- 
element 
solution 

0.5 

0.4438 

0.4390 

- 0.07137 

0.6 

0.4419 

- 0.06208 

■ 

0.7 

0.4419 

- 0.05165 

- 0.05302 

0.005925 

0.8 

0.4413 

0.4413 

- 0.04489 

0.9 

0.4404 

- 0.03964 

- 0.04145 

0.01303 

1.0 

0.4389 

0.4392 

- 0.03546 

- 0.03712 

0.01780 

0.01761 

1.2 

0.4362 

0.4360 

- 0.02925 

- 0.03060 

0.02838 

1 .4 

0.4332 

0.4334 

- 0.02487 

- 0.02592 

0.03961 

0.03917 

1.6 

0.4302 

- 0.02163 

1 .8 

0.4265 

- 0.01917 

- 0.02029 

0.06032 

0.05967 

2.0 

0.4228 

0.4229 

- 0.01815 

- 0.01818 

0.06894 

0.06826 

0.5 

0.3650 

0.3719 

- 0.1310 

0.001902 

0.001886 

0.6 

0.3650 

0.3653 

E 

- 0.1277 

0.003672 

0.003648 

0.7 

0.3638 

0.3642 

- 0.1139 

- 0.1171 

0.006227 

0.006175 

0.8 

0.3638 

0.3632 

- 0.1060 

- 0.1085 

0.009542 

0.009448 

0.9 

0.3622 

0.3626 

- 0.1003 

- 0.1039 

0.01348 

0.01333 

1.0 

0.3622 

0.3618 

- 0.09605 

- 0.09925 

0.01783 

0.01762 

1.2 

0.3594 

0.3603 

- 0.09029 

- 0.09415 

0.02680 

0.02646 

1.4 

0.3573 

0.3583 

- 0.08670 

- 0.08958 

0.03497 

0.03451 

1 .6 

0.3550 

0.3560 

- 0.08432 

- 0.08628 

0.04167 

0.04112 

1.8 

0.3525 

0.3541 

- 0.08268 

- 0.08070 

0.046*0 

0.04625 

2.0 

0.3498 

0.3541 

- 0.08150 

- 0.07757 

0.05090 

0.05021 
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conditions  (i.e.,  simply-supported),  if  the  loading  is  changed  to  that  of 
uniform  loading,  closed-form  solution  does  not  exist.  However,  the  finite- 
element  method  can  predict  the  neutral-surface  locations  (which  may  now 
depend  on  x and  y)  and  deflections. 

Table  4.5  shows  the  nondimensional  neutral -surface  locations  at  the 
center  of  the  plate  (i.e.,  in  element  1),  and  deflections  for  the  problem 
of  Table  4.4.,  but  with  uniform  loading.  The  trend  and  magnitudes  of  the 
solutions  are  similar  to  those  in  Table  4.4  for  the  sinusoidal  loading. 
Table  4.6  shows  the  nondimensional  neutral-surface  locations  for  various 
elements.  It  is  clear  that  zx  does  not  vary  much  along  x or  y.  However, 
z varies  noticeably  with  y while  it  is  almost  constant  with  respect  to  x. 
Also,  the  locations  change  with  the  plate  aspect  ratio. 

As  pointed  out  earlier,  the  present  finite- element  formulation  is  also 
good  for  thick  plates  (i.e.,  accounts  for  transverse  shear  deformation). 
Representative  results  are  presented  to  show  the  effect  of  the  thickness 
on  the  neutral-surface  locations  and  deflections  (see  Reddy  and  Bert  [15]). 

Figure  4.1  shows  the  influence  of  the  aspect  ratio  (b/a)  and  side-to- 
thickness  ratio  (a/h)  on  the  location  of  neutral  surfaces  for  a single-lay- 
er, orthotropic,  bimodulus,  simply  supported  rectangular  plate  subjected 
to  sinusoidal  loading.  The  following  elastic  properties  are  used: 

Eli  = 3.584  GPa  , E^  = 1.792  GPa  , e£2  = Eu  , E§2  = E^ 

^12  = G?2  = 1-27  GPa  , vj2  = v21  = 0.4  , vj2  = v21  = 0.2 


Note  that  for  b/a  = 1,  the  neutral  surfaces  associated  with  x-  and  y- 

directions  coincide  (i.e.,  z_  = z ). 

' nx  ny 
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Table  4.5.  Neutral -surface  locations  and  deflections  for  simply- 
supported  rectangular  plates  of  cross-ply  laminated 
aramid-rubber  and  polyester-rubber  (uniform-loading  case) 


0.9  : 0.3622 

1.0  ; 0.3612 

1.2  0.3590 

1.4  0.3564 

0.3533 
0.3493 
0.3445 


0.19040 

0.14270 

0.12160 

0.10980 

0.10230 

0.09714 

0.09046 

0.08648 

0.08390 

0.08182 

0.08043 


j 0.00287 
0.00568 
0.00976 
I 0.01507 
; 0.02137 
| 0.02832 
0.04258 
0.05544 
0.06580 
| 0.07361 
I 0.07927 


2.3433 

1.9502 

1.6709 

1.4622 

1.3000 

1.1724 

0.9818 

0.8495 

0.7543 

0.6844 

0.6332 


0.5546 
0.4555  I 
0.3960 
0.3583  i 
0.3330 
0.3152 
0.2918 
0.2782 
0.2693  ' 
0.2624  i 
0.2571  ! 
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Table  4.6.  Neutral -surface  locations  as  a function  of  the  element 
location  for  simply-supported  rectangular  plates  of 
cross-ply  laminated  aramid-rubber  and  polyester- 
rubber  (uni form- loading  case) 


Mat- 

erial 


Element 

No. 


a/b  = 0.5 


0.44397 

0.44398 

0.44401 

0.44405 

0.44407 

0.44325 

0.44197 

0.43980 

0.43796 


- 0.11414 

- 0.11481 

- 0.11618 

- 0.11792 

- 0.11919 

- 0.08752 

- 0.05971 

- 0.04162 

- 0.03378 


0.44002 

0.44013 

0.44033 

0.44059 

0.44079 

0.43947 

0.43835 

0.43686 

0.43572 


0.043286 

0.043781 

0.044708 

0.045820 

0.046613 

0.040654 

0.036348 

0.031648 

0.028602 


0.36597  1-0.19040 

0.36598  1-0.19101 


0.36600 

0.36603 

0.36605 

0.36526 

0.36415 

0.36290 

0.36199 


- 0.19224 

- 0.19386 

- 0.19505 

- 0.15496 

- 0.12374 

- 0.10556 

- 0.098157 


0.36124 

0.36137 

0.36161 

0.36190 

0.36211 

0.36101 

0.36055 

0.35999 

0.35961 


0.097142 

0.097651 

0.098586 

0.099683 

0.010045 

0.096155 

0.094110 

0.091909 

0.090318 


0.41739 

0.42006 

0.42421 

0.42841 

0.43094 

0.41693  j 

0.41612 

0.4523 

0.41463  i 

i 


0.34451  ! 

0.34755  ■ 

0.35165 

0.35519 

0.35709 

0.34431 

0.34397 

0.34358 

0.34331 


0.016405 
0.017267 
0.018940 
0.021052 
0.022623 
0.016184  I 
0.015769  ! 
0.015219  ; 

0.014782  I 

■ i 

i 


0.080437 
0.081660 
0.082619  i 
0.084542  j 
0.086014  | 
0.080356  j 
0.080201  j 
0.079991 
0.079821 


18 


Similar  results  are  presented  in  Figs.  4.2  and  4.3  for  a two-layer, 
cross-ply  (0°/90°),  rectangular  plate  of  aramid-rubber  under  sinusoidal 
loading.  Note  from  Fig.  4.2  that  the  neutral -surface  location,  Zx  is 
virtually  unchanged  for  aspect  ratios  greater  than  1,  while  the  neutral- 
surface  location,  Z^,  increases  proportional  to  the  aspect  ratio.  It 
should  also  be  noted  that  the  neutral  surfaces  do  not  coincide  in  the  cross 
ply  case  for  b/a  = 1 . 

Figure  4.4  shows  the  influence  of  the  aspect  ratio,  and  side-to-thick- 
ness  ratio  on  the  transverse  deflection  for  the  single-layer  and  two-layer 
cross-ply  plates  discussed  above.  The  effect  of  thickness  on  the  deflec- 
tion is  more  pronounced  than  the  effect  of  the  aspect  ratio. 

5.  SUMMARY  AND  CONCLUSIONS 

A penalty  finite-element  formulation  of  the  equations  governing  the 
small  deflection  elastic  behavior  of  thin  plates  laminated  of  anisotropic 
bimodulus  materials  is  presented.  The  resulting  finite  element  model  is 
valid  for  the  analysis  of  thin  and  moderately  thick  plates.  Single-layer 
and  two-layer  cross-ply,  simply-supported  thin  rectangular  plates  subjected 
to  sinusoidally  distributed  normal  pressure  and  uniformly  distributed  normal 
pressure  are  analyzed.  In  the  case  of  sinusoidal  loading,  the  present 
finite-element  solutions  agree  very  closely  with  the  exact  closed-form  solu- 
tions. To  illustrate  the  validity  of  the  present  finite  element  model  for  thick  plates 
numericalresults  are  presented  for  side-to- thickness  ratio  of  10  for  single- 
layer and  two-layer  cross-ply  rectangular  plates. 

Acknowledgement.  The  results  reported  in  this  paper  were  obtained  during 
an  investigation  supported  by  the  Office  of  Naval  Research,  Structural 
Mechanics  Program  through  Contract  N00014-78-C-0647.  The  authors  wish  to 
express  their  sincere  thanks  to  Professor  C.  W.  Bert  for  helpful  discussions. 
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